Oscillations of Andreev states in clean ferromagnetic films 
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We investigate the influence of the exchange field on the 
Andreev bound states in a ferromagnetic (F) film backed on 
one side by a superconductor (S). Our model accounts for 
diffusive reflection at the outer surface and possible backscat- 
tering at the FS-interface. Phase shifting of the Andreev level 
by the exchange field results in an oscillatory behaviour of the 
density of states of F as a function of the layer thickness. We 
show that our results agree quantitatively with recent exper- 
iments. 

Probing the proximity effect by tunneling spectroscopy 
of induced superconducting correlations has a long his- 
tory. n Early experiments on the proximity density of 
statesEl could be understood in the tunneling model of 
McMillanB. Recently the spatial dependence of the 
proximity- density of states in normal metals has been 
measuredQ and successfully explained in terms of the qua- 
siclassical theoryQ. The influence of spin-splitting by a 
parallel magnetic field was measured in ^| and found to 
coincide with a Zeeman-split density of states. A new 
experimental and theoretical challenge is to extend these 
studies to the superconducting proximity effect in ferro- 
magnets. 

Previous experimental investigation have concentrated 
on thermodynamic properties of FS- multilayers. Here os- 
cillations of the superconducting critical temperature T c 
as a function of the thickness of the FJayers have been 
predicted!! and studied experimentallyQ. However, the 
experimental evidence for these T c -oscillation is not con- 
clusive. This may result from pan incomplete knowledge 
about the FS-interface qualityu In our opinion, it is also 
questionable if the theoretical approach using the diffu- 
sive quasiclassical formalism^ is applicable for F-layers of 
typical thicknesses 10 — 50A. 

The most recent experiments have concentrated on 
other properties of FS-layers. Ryazanov et altB studied 
the supercurrent through a thin ferromagnetic layer and 
found a non-monotonic temperature dependence, which 
can be interpreted in terms of a -«j>hase shift due to the 
exchange splitting. Kontos et alt3, on the other hand, 
have studied the density of states in a thin ferromagnetic 
layer in contact to a superconductor. An oscillatory be- 
haviour of the induced superconducting correlation was 
observed for layers of different thicknesses and attributed 
to influence of the exchange field. It is this experiment, 
that motivated our present study. 

Bearing this in mind, we study the superconducting 
proximity effect in a thin ferromagnetic layer. The F-film 



is characterized by an homogenous exchange splitting h. 
We model the film as a ballistic layer with rough bound- 
aries. Band mismatch and disorder at the interface may 
lead to enhanced backscattering at the FS-boundary. We 
will derive a general formula for the subgap density of 
states depending only on the length distribution of clas- 
sical trajectories in the F-layer. The resulting density 
of states shows as a signature of the exchange splitting 
an oscillatory behaviour as a function of layer thickness. 
Comparison to the experimental data show reasonable 
agreement keeping in mind the large uncertainty of some 
sample parameters. 




FIG. 1. a) Schematic drawing of a ferromagnetic film in 
connection with a superconductor. A typical classical trajec- 
tory is also indicated. An electron coming from the bulk of 
the superconductor enters into the ferromagnet (at point i) 
and after several diffusive reflections from the insulator and 
the SF- interface returns to the superconductor (at point o). 
b) The calculated distribution of the trajectory lengths in the 
F-layer for small SF-transparency (here T — 0.1). The dou- 
ble peak structure close to the smallest length originates from 
the first two reflections from the insulator, whereas the dis- 
tribution for long trajectories decays as exp(— l/l), with the 
average length I ~ 2d/T 

The system we consider is sketched in Fig. [l| A thin 
ferromagnetic layer (F) of thickness d is connected to a 
superconducting bank (S) on one side and bounded on 
the other side by an insulator. The F-layer is charac- 
terized by an exchange splitting, which we take into ac- 
count as mean field h in the Hamiltonian. The thickness 



1 



d is larger than the Fermi wave length Af and smaller 
than the elastic mean frae-,path £- lmp which allows for a 
quasi-classical description^ in the clean limit. Then, the 
real-time Eilenberger equation reads 

-iv F Vg a {E,v F ,r) = (1) 
[(E + ah(r))r 3 - if 2 A(r), g a (E, v F , r)} . 

Here f< denote the Pauli matrices, A(r) is the (real) su- 
perconducting pair potential, and a (= ±1) labels the 
electron spin. The matrix Green's functions have to obey 
the normalization condition g%. = 1. Inside the F-layer 
h is constant and A = 0. We neglect the change of the 
pair potential in S, thus A(r) = const, inside the super- 
conductor. We have solved Eq. (|l|) along each classical 
trajectory in F that comes from the superconductor and 
ends there. The density of states for \E\ < A on a given 
trajectory of length I is then given by 

N(E,l) = ^ ^e[TrT 3 g a (E,v F ,r)} 



<T=±l 
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Here No is the density of states at the Fermi level in 
the normal state. Inside the F-layer N a (E, I) is constant 
along a given trajectory and depends only on the length 
of the trajectory I. Eq. (||) means that the density of 
states below A is a sum of (5-peaks resulting from Andreev 
bound states. The energies follow from the quasi-classical 
quantization I — l n . The total density of states can then 
be found by averaging (||) over all classical trajectories. 
Denoting the trajectory length distribution p(l) we find 
for the density of states 
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This formula presents the general result for subgap den- 
sity of states of a quasiballistic ferromagnet connected to 
a superconductor. It is completely specified by the length 
distribution of classical trajectories. This depends only 
on the geometrical properties of the attached ferromag- 
net and the connecting interface. 

Now we have to specify the trajectory length distribu- 
tion for our particular case. We model the thin layer by 
a weakly disordered thin film with a rough surface and 
a rough SF-interface of average transparency T. A typi- 
cal trajectory is depicted in Fig. 1. An electron coming 
from the bulk of S enters into the F-layer and after sev- 
eral reflections from the insulator and the FS-interface 
returns to the S-bank (see Fig. |l|), where it is Andreev 
reflected as a hole which traverses the trajectory in the 



opposite direction. Thus, the elementary building block 
of a typical trajectory is the segment between two suc- 
cessive reflections from the superconductor. The number 
of blocks which form the total trajectory depends on the 
transparency of the interface, i. e. it is roughly ~ 1/T. 
First, let us consider the length distribution of one el- 
ementary block. Due to the roughness of the insulator 
and the FS-interface the quasi-particles undergo diffusive 
reflection from these boundaries. Incoming and outgoing 
direction are completely uncorrelated. Accounting for 
the weak bulk disorder we include a factor exp(— l/£- lmp ). 
The length distribution of one elementary block is then 
given by 
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where C = E^d/^imp) (E 2 (z) = dxexp (-zx)jx 2 is 
the exponential integral of second order). Second, we 
connect the elementary building blocks, if the S-F inter- 
face has an average transparency T. In determining the 
length distribution we assume that an particle either goes 
through the interface or is fully reflected. Only the num- 
ber of these reflections depend on T. We do not take into 
account quantum mechanical interference for a single re- 
flection at the FS-interface. Taking this into account will 
yield essentially the same results as our approach. By an 
expansion in the reflectivity R = 1 — T it is easy to see 
that the full distribution p(l) obeys the integral equation 



p(l)=T Po (l)+R / dl'p (l') P {l-l') 



(G) 



This is readily solved by a Fourier transformation. We 
obtain 
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where Po(k) = Fi^iikd + d/£i mp )/C is the Fourier trans- 
form of Po(l). The distribution p(l) is plotted for T = 0.1 
and d/£ imp — 0.1 in Fig. [l| It's characteristics are a 
double peak structure close to the shortest trajectories 
I > 2d resulting from trajectories reflected once or twice 
from the insulator. At large I the distribution decays ex- 
ponentially as exp(— l/l), where I = 2d\n(£ lmp / d) /T is 
the mean trajectory length. For T ~ 1 p(l) has only one 
peak close to 2d. We therefore have two characteristic 
lengths of the distribution, the smallest possible trajec- 
tory length 2d and the average length I. Which of these 
length scales will determine the total density of states 
will depend on the other parameters, in particular on h. 

Combining Eqs. (§, (§ and (@) we obtain for the total 
DOS 
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where k n — 2n(E + ah)/vp. Thus, the density of 
states is fully expressed in terms of known functions. 
In the following we will discuss the parameter range of 
rather strong exchange fields in the limit of thin layers 
d -c vp/A. The results presented below are in the most 
realistic case of small SF- interface transparency T< 1. 
In this limit the distribution of long lengths is well ap- 
proximated as exp(— I /I). Most probably the F-film has 
a non-uniform thickness due to the large scale roughness 
of the boundaries. For a smoothly varying thickness we 
can take this into account by averaging the (^) over a 
Gaussian distribution of the thicknesses around a mean 
value d. This will also lead to a smoothening of the sharp 
features in the DOS resulting from the lower cutoff in 
p{l). The qualitative behavior will however not change. 
In practice, we have chosen a width of the distribution 
to be of order 10% that corresponds to condition of the 
experiment!!^. 



(nvp/hl) exp(— 7TOf/2W)/(1 — exp(— Ttvp/hl)), which fol- 
lows from the approximation mentioned above. At finite 
energies N(E) passes through minima corresponding to 
E = ±h before approaching the normal metal DOS. 

For larger exchange fields when hd > 1 mainly the 
short trajectories of I ~ 2c? contribute to the energy de- 
pendence of N(E). The DOS is close to that of the nor- 
mal state. The interesting part are the small deviations 
of amplitude T which oscillate as a function of dh. In 
Fig. 3 we have plotted (N(E)/N - 1)/T versus dh/v F 
for different energies. At small hd/vp <C 1 and E -C A 
the DOS as a function oihd/vp still reminds of the double 
peak structure of the length distribution. For hd/vp ^ 1, 
however, it develops coherent oscillations as a function of 
hd/vp with a period of 7r/2. Note that the magnitude and 
the sign of the oscillation depend on the energy. Maximal 
amplitudes of opposite sign always correspond to E = 
or \E\ = A. 
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FIG. 2. Density of states for weak exchange fields and 
small SF-interface transparency. For E = the curve is in- 
dependent of T as long as T <C 1. At finite energies the 
precise form depends on the transparency (here T = 0.03 and 
^im P = lOOd). It has a broad peak at zero energy and two 
dips E — ±h resulting from the shifts of the Andreev bound 
states. 

For weak exchange fields when hl/vp < 1 the DOS 
at low energies E <C A is mainly contributed from long 
trajectories with I > /. The DOS versus hd/T is plot- 
ted for different energies in Fig 2. For h = (a normal 
film), N(E) vanishes at zero energy and increases for fi- 
nite energies ending with a peak at A. The effect of the 
exchange field in the regime h < A is to split the DOS 
for the two spin bands. Thus the total DOS is the aver- 
age of two by ±/i shifted spectra resulting in a peak at 
zero energy and two dips at E = ±h. Increasing h fur- 
ther leads to a suppression of superconducting features 
of the DOS. The zero energy DOS increases roughly as 
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FIG. 3. Oscillation of reduced DOS (N(E)/N - 1)/T 
with the exchange field. For hd/vp > 1, it develops coherent 
oscillations as a function of dh/vp with a period of n/2 and an 
energy dependent magnitude and sign. Maximal amplitudes 
of opposite sign always correspond to E = or \E\ — A. 

Let us now discuss the relation of our results to the 
recent experiment performed in the group of AprihO. 
They observed DOS oscillation in thin films of the ferro- 
magnetic alloy Pdi_ x Ni x with x of the order of 10% in 
contact with a Nb electrode. By tunneling spectroscopy 
they measured the differential conductance versus bias 
of the F-film at low temperatures ~ 300mK. Normalized 
to the normal state conductance this yields the density 
of states in the ferromagnetic layer. For two different 
thicknesses of the F-film, d\ — 50A and g?2 = 75A the 
reduced density of states differs in sign and magnitude, 
but has the same functional energy dependence. The 
exchange field was estimated by several methods to be 
in the range h = 5 — 15meV. These experimental deter- 
minations were, however, only sensitive to the average 
magnetization. In thin layers, it is reasonable to assume 
that the magnetic structure has multiple domains, inho- 
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mogeneous thickness (as already discussed before), and 
non-uniform Ni concentration. All these effects tend to 
suppress the measured values of the exchange field below 
the value, which results from the local Ni-doping level. 
We will therefore regard the experimental values of h as 
lower bounds to the value used in our fits. In Fig. 4 we 
compare the data of RefO with our calculation. We have 
taken v F = 2x 10 7 cm/s and h = 29meV, and T = 0.055. 
The agreement between the experimental data and our 
calculation is satisfactory. Both, the functional depen- 
dence and the sign change are correctly reproduced by 
our calculations. 




FIG. 4. Energy dependence of DOS for two different thick- 
nesses of the ferromagnetic film, di — 50A (solid line) and 
c?2 = 75 A (dashed line), with the corresponding experimental 
curves of RefJid. The DOS oscillation appears as an inversion 
of the energy dependence of DOS when d changes by 7tvf/4/i. 
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In conclusion we have investigated theoretically the su- 
perconducting proximity effect in thin ferromagnetic lay- 
ers in a quasi-ballistic model. We have found that the 
effect of the ferromagnet exchange field h is to suppress 
the superconducting features in the density of states for 
dh/v F > T. At exchange fields larger than qvp/d the 
density of states oscillates around the normal state value 
as a function of dh/v F with a period of 7r/2 and an am- 
plitude of the order of the interface transparency T. We 
have shown that the oscillation can lead to an inverted 
energy dependence of the density of states.— This effect 
has been observed in the experiment of RefO. We have 
achieved quantitative agreement between our theory and 
the experimental data. 
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